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ABSTRACT 

We calculate the degree of linear polarization of radiation from stars having 
planets that may not be resolved spatially. We assume single scattering by water 
and silicate particulates in the planetary atmosphere. The dilution of the re- 
flected polarized radiation of the planet by the unpolarized stellar radiation and 
the effect of oblateness of the planet as well as its elliptical orbit are included. 
We employ a chemical equilibrium model to estimate the number density of wa- 
ter and silicate condensates and calculate the degree of linear polarization at R 
band of starlight as a function of (1) mean size of condensates, (2) planetary 
oblateness, (3) inclination angle, (4) phase angle, (5) orbital eccentricity e and 
(6) the epoch of periastron passage. We show that the polarization profile alters 
significantly at all inclination angles when elliptical orbit is considered and the 
degree of polarization peaks at the epoch of periastron passage. We predict that 
detectable amount of linear polarization may arise if the planetary atmosphere is 
optically thin, the mean size of the condensates is not greater than a few microns 
and the oblateness of the planet is as high as that of Jupiter. 

Subject headings: binaries: general - dust, extinction - polarization - scattering - 
planetary systems 
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1. INTRODUCTION 

Polarization has always been an efficient tool to probe the physical properties in the 
environment of various astrophysical objects. Recently, it is realized that polarization could 
be a very important diagnostic method for analyzing the atmosphere of sub-stellar massive 
object such as brown dwarfs and extra-solar planets. Sengupta & Krishan (2001) predicted 
detectable amount of linear polarization from L dwarfs because of the presence of condensates 
in the visible atmosphere and subsequently linear polarization at optical is detected by 
Menard. Delfosse & Monin (2002) and Zapatero Osorio et al. (2005). The observed linear 
polarization can well be explained by single dust scattering models that assumes optically 
thin and rotationally oblate photosphere (Sengupta 2003; Sengupta & Kwok 2005). 

On the other hand, the use of polarimetry in detecting and understanding the physical 

properties of extra-solar planets is emphasized by Seager et al. (2000); Saar & Seager (2003); 
Stam et al. (2004). While Seager et al. (2000) modeled the polarization due to close-in- 
planet or the so called roaster such as the fcst discovered extra-solar planet 51 Peg b (Mayor 
& Queloz 1995), Stam et. al (2004) estimated the degree of polarization caused by planet 
at Jupiter's distance from the Sun. In both the methods, the reflected flux of the planet 
and the Stokes vectors of the reflected radiation are calculated by assuming circular orbit 
and blackbody radiation from the star. The planet is considered to be spherical and the 
polarization is integrated over the illuminated portion of the disk. While multiple scattering 
polarization is calculated by Stam et al. (2004), Seager et al. (2000) employed Monte-Carlo 
proceedure and computed the number of scatterings. The later authors found that for 
an almost absorbing atmosphere, higher amount of polarization may arise due to single 
scattering. In both the investigations, the observable degree of polarization is estimated by 
simply multiplying the flux ratio with the polarization of the radiation from the spatially 
resolved planet. Since degree of polarization is a relative measure, it does not carry any 
information on the radius of the circular orbit or the planetary radius. 

However, in the most realistic situation, the planet should be oblate because of rotation 
around its own axis. The oblateness of solar planets ranges from 0.003 for Earth to 0.065 
for Jupiter and 0.1 for Saturn. Further, most of the extra-solar planets detected so far have 
eccentric orbit around the star and the eccentricity ranges from 0.0 to as high as 0.7. 

Although multiple scattering is a reasonable choice for optically thick medium, it un- 
derestimates the amount of polarization by a few orders if the medium is optically thin and 
hence if polarization is caused by single scattering. As mentioned before, single scattering 
polarization could successfully explain the high degree of linear polarization observed from L 
dwarfs. Therefore, it's quite reasonable to investigate the amount of polarization that could 
arise due to single scattering in extra-solar planets. 



In this paper, we present polarization profiles of starlight caused by single scattering in 
the atmosphere of an oblate planet rotating in elliptical orbit and show that in certain cases 
the polarization may be detectable by the existing instrumental facilities. 



Since the dust density is assumed to be low and scatterings by atoms and molecules 
(e.g., Rayleigh scattering) do not contribute to polarization significantly, single scattering 
approximation is reasonable for the region where the optical depth r < 1. If present, multiple 
scattering can reduce the degree of polarization by a few orders of magnitude (Sengupta 
& Krishan 2001) because the planes of the scattering events are randomly oriented and 
average each other's contribution out from the final polarization. Hence, the amount of 
observed linear polarization could act as a probe in favour of single or multiple scattering 
approximation. 

In the present work, we use the formalism given by Simmons (1983) which is a general- 
ization of the work by Brown, McLean & Emslie (1978). In this formalism, the primary star is 
assumed to be a point source of unpolarized light. For most practical application, the Stokes 
parameter are normahzed and defined by [Irefi, Q, U, V]/Itotai where Itotai = ^refi + Istar, hefi 
being the reflected flux or intensity from the planet and Igtar being the unpolarized flux or 
intensity received from the star. We have neglected the thermal radiation of the planet if any 
because for a sufficiently old planet the thermal radiation should be much less as compared 
to the reflected radiation or there may not be thermal radiation by comparatively smaller, 
earth like planets. However, for young giant planets far away from the primary star, thermal 
radiation should be significant. Since L,.f.fi is much less than I star-, "wc consider Itotai = I star 
and hence after we shall refer Irefh Qi U and V to normalized Stokes parameter. Hence the 
degree of linear polarization can be written as P = ^/Q^^^^[P 

The density distribution of scattcrcr is calculated in co-rotating reference frame so that 
it is independent of the phase angle A of the planet. The transformation to observer frame 
is done by using the properties of spherical harmonics under rotation. In a circular orbit, 
the normalized Stokes parameter Q and U are given as a harmonic series : 



2. 



SINGLE SCATTERING POLARIZATION OF STARLIGHT 
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where k — 27r/A, A being the wavelength, i is the orbital inclination angle. The harmonic 
co-efRcients are given by 
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M — max{2,m) and Hl^{i) are given in Simmons (1983). rjim and are related 

with the density distribution in the co-rotating frame and are given by 
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a(l,m) f n' (r,e,(l))Pr( cos 0i) ( ) sinddedctxir. (5) 
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where n'(r, 6,(f)) is the number density of scatterer in the co- rotating frame, 9i is the viewing 
angle, PJ^ is the associated Legendre function of the first kind and 
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At an edge-on view, 9i = 7i/2 and = 0, and hence $,ijn — 0. Fi2{k) is related to the 
scattering function and is given by 



Im 



a{l,m) f ''^^''^^ ''^^''•^^ P;"'(cosg)d(cosg). (7) 
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In the above equation, 6 is the scattering angle, ii and 12 are the scattering functions given 
by van de Hulst (1957). ii and 12 depend on the refractive index as well as on the size and 
shape of the scatterer. In the present work we consider spherical dust particles as scatterer. 

Considering an ellipsoidal density distribution and using the addition theorem of spher- 
ical harmonic, rjim can be written as 



?7i^ = 27ra(Z,m)Pr(0) r\{r)dr I 

JR2 J-i 



where Ri and R2 are the outer and the inner equatorial axis length of the planet, A is the 
ratio of the length of the equatorial axis to the polar axis, and n = cos 6. We have taken 
multi-poles up to 1=5. We have considered up to fifth harmonic, i.e., m=0,l,2,3,4,5. However, 
since the density distribution is symmetric about the orbital plane, the first, the third and 
the fifth harmonics are zero and the fourth harmonic is small as compared to the second 
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harmonic. In other word, for the adopted density distribution, the odd number harmonics 
are zero and the degree of polarization is determined mainly by the second harmonic. 

We convert n{r)dr into n{P)dP by using the equation of hydrostatic equilibrium n{r)dr — 
n{P)dP/ gp{P), where P is the pressure at different geometrical height, p is the mass density 
at different pressure scale, and g is the surface gravity (which can be assumed to be constant 
for a geometrically thin atmosphere). The degree of polarization is calculated at wavelengths 
ranging from 0.5486 to 0.8491 //m and averaged by using the response function of R band 
Bessel filter. 



3. EFFECT OF ORBITAL ECCENTRICITY 

The effect of orbital eccentricity enters through a multiplicative factor /i(A, e) to the 
co-efficient of the dominant harmonic p2 and V2 (Brown et al. 1982) and it is given by 

^^^^^^Al±I0>_z^ (9) 

where e is the orbital eccentricity and A is the true anomaly which is related to the eccentric 
anomaly E by the relationship 

tan^^j^^^j tan-, (10) 

and Ap is the longitude of the periastron. The eccentric anomaly E is related to the orbital 
phase angle through Kepler's equation 

£;-esin£; = A- Ap (11) 

where A — Ap = 27r(t — Tq)/P is the mean anomaly, t being any epoch of time, P is the 
orbital period and Tq is the epoch of periastron passage. We have assumed no drastic seasonal 
change in the planetary atmosphere and that the atmospheric T-P profile remains the same 
through out the planetary year. 

For a circular orbit, the polarization profile is determined by the second harmonic only. 
But for an elliptical orbit, the first and the third harmonics are non-zero although the second 
harmonic remains dominant. 



4. THE ATMOSPHERIC MODELS 



We have adopted the atmospheric models of extra-solar planets given in Sudarsky et al. 
(2003). The Temperature-Pressure (T-P) profiles are taken for models of a class II and a 
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class V planets or closed-in-planet also known as roaster. The atmosphere of class II planets 
should have water as condensates while class V planets should contain silicate in their upper 
atmosphere. In calculating the T-P profile of the planet Ups And d (a class II planet), 
Sudarsky et al. (2003) considered the surface gravity to be 2 x 10^ cms^^ while for the 
planet HD 209458 b (a class V planet) the surface gravity is taken to be 980 cms~^. In the 
present work, we have taken the same values for the surface gravity of the two representative 
planets. We have assumed that the thermal radiation of the planet is negligible compared to 
the reflected radiation so that the contribution to polarization comes only from the reflected 
radiation. 

The dust distribution in the atmosphere is calculated based on the one dimensional 
cloud model of Cooper et al. (2003). This model assumes chemical equilibrium throughout 
the atmosphere, and uniform density distribution across the surface of an object at each 
given pressure and temperature. Under these assumptions, the number density of cloud 
particles is given by 



where p is the mass density of the surrounding gas, r is the cloud particle radius, is the 
mass density of the dust condensates, p and p^ are the mean molecular weight of atmospheric 
gas and condensates respectively. The condensate mixing number ratio {qc) is given as qc — 
QbeiowPc,i/P for heterogeneously condensing clouds where qbeiow is the fraction of condensible 
vapor just below the cloud base, Pc,i is the pressure at the condensation point, and P is the 
gas pressure in the atmosphere. The condensation curves for water and forstcritc condensates 
arc taken from Sudarsky ct al. (2003). The values of pd, Pd and qheiow for forsterite and water 
are taken from Cooper ct al. (2003). The real part of the refractive index for forsterite is 
fixed at 1.65 at any wavelength and the imaginary part is taken by interpolating the data 
given in Scott & Duley (1996). For the refractive index of water, we have used the data 
given in Segelstein (1981). 

Apart from the calculation of the grain number density, the location of the cloud in the 
atmosphere plays an important role in determining the amount of polarization. The location 
of the cloud base for different atmospheric models and different chemical species is determined 
by the intersection of the T-P profile of the atmosphere model and the condensation curve 
Pc,i as prescribed in Cooper et al. (2003). Taking the condensation curve for forsterite and 
water as given in Sudarsky et al. (2003), we determine the base of the cloud from the T-P 
profiles of Ups And d and HD209458 b calculated by the same authors. 

At present, there is no convincing justification in favor of any specific form of the particle 
size distribution function. In the present work, we adopt a log-normal size distribution 
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function used by Ackerman & Marley (2001) and Saumon et al. (2000) with a fixed width 
at 1.3. Following Sudarsky et al. (2003), we have taken the mean particle radius for water 
5.0 /im and that for silicate (forsterite) 10.0 /xm. However, in order to show the effect of 
mean grain size on the degree of polarization we have also presented polarization profile 
as a function of mean particle size keeping the T-P profile the same as that calculated by 
Sudarsky et al. (2003) for water with mean grain size 5.0 ^m. 

In multiple scattering, the mean particle size should be different at different pressure 
scale. However, for single scattering, a fixed mean particle size of a particular condensate is 
sufficient (Sengupta & Kwok 2005). 



The rotationally induced oblateness of solar planets has been discussed in details by 
Hubbard (1984) and Murray & Dermott (2000). Recently, the formalism for oblateness 
is extended to extra-solar planets by Barnes & Fortney (2003) who used Darwin-Radau 
relationship 



(Murray & Dermott 2000; Barnes & Fortney 2003) to relate rotation to oblateness. In the 
above equation, K = I/MRl < 2/3 is the moment of inertia parameter of an object with 
moment of inertia /. The Darwin-Radau relationship is exact for uniform density objects 
{K — 0.4) and provides a reasonable (within a few percent of errors) estimation of the 
oblateness of the solar planets. At 1 bar pressure level, the oblateness / of Jupiter, Saturn and 
Earth are 0.065, 0.098 and 0.003 respectively. Barnes & Fortney (2003) modeled the planet 
HD209458b and estimated its oblateness to be about 0.00285. Apart from rotational effects, 
tidal interaction of a closed in planet with its primary star may also impose an ellipsoidal 
shape extending toward the star. In the present work we ignore such effects. Moreover, the 
estimation of the moment of inertia of the planets needs density distribution, radius and 
mass of the planet along with its rotational period and so it is highly model dependent. In 
the present work we adopt a wide range of values for oblateness, the maximum value is equal 
to that of Jupiter (/ = 0.065) and the minimum value is equal to that of Earth(/ = 0.003). 



5. 



THE OBLATENESS OF PLANETS 
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6. RESULTS AND DISCUSSION 

We present in figure 1 and in figure 2, the polarization profiles at R band (Bessel filter) 
of star with a planet having forsterite as a major condensates in its atmosphere. Silicate and 
iron can form high in the atmosphere of class V roasters that orbit their stars around 0.05 
AU. Although presence of more than one species affects the amount of polarization because 
of differences in size, optical properties, and location, we assume forsterite as the major 
condensate in our model for simplicity. However, the change in the amount of polarization 
by the incorporation of other condensate species cannot be determined due to the domi- 
nant effect of other parameters such as inclination angle, orbital eccentricity, oblateness etc. 
Figure 3 presents the degree of polarization as a function of planetary oblateness. Figure 4 
presents the polarization profile for planets with non-circular orbit and the change in degree 
of polarization caused by the change in orbital eccentricity is shown in figure 5 . 

Our polarization profiles for the orbital eccentricity e=0 i.e., for circular orbit, presented 
in figure 1 and figure 2, are qualitatively the same as that presented by Stam et al. (2004). 
When the inclination angle i = 0, the observed part of the planet that is illuminated is con- 
stant and hence the degree of polarization is constant. The maximum amount of polarization 
produced by the planet is the same for all values of i. When the inchnation angle i — 90°, 
the polarization is zero at the fractional orbital period t/P — 0.25 and 0.75 that correspond 
to the phase angle A = 180° and 0° respectively because the planet's night and day side are 
turned toward the observer in these cases. Apart from the geometrical asymmetry of any 
object, another asymmetry plays important role for a binary star or a planet around a star 
in determining the polarization profile. This second effect is due to the asymmetry in the 
position of the object (or the scatterer) with respect to the observer and the source of light. 
As the symmetry axis of the planet changes position with respect to the observer, the total 
number of effective scatterers also changes. This effect can be visualized as the amount of 
light that is scattered towards the observer from the net illuminated region. In a circular 
orbit, the observed degree of polarization is determined by the scatterers that are present 
within the illuminated area facing towards the observer. The observed degree of polarization 
is zero when the night side of the object is towards the observer and it increases as the 
illuminated area towards the observer increases. The fraction of the area that is illuminated 
depends on the inclination angle as well. As the inclination angle decreases, the orbital plane 
becomes more and more face on towards the observer. Consequently, the illuminated area 
towards the observer increases yielding into higher amount of polarization that can be seen 
by the observer. 

However, the polarization profile alters drastically when the orbit is non-circular as can 
be seen in figure 4. The degree of polarization peaks when the planet passes its periastron. 
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Hence, the degree of polarization is time variable even if the inclination angle is zero degree. 
The change in the amount of polarization is determined by the eccentricity of the orbit as 
can be seen from Figure 5. 

As discussed by Brown et al (1982), the "weighted" or effective optical depth (as seen 
by the source) of the scattering region idealized as point scatterer containing N particles is 
proprtional to the square of the ratio between the semi-major axis a and the stellar separation 
at longitude A given by -R(A) = a(l — e^)/[l + e cos(A — Ap)]. For a circular orbit this ratio is 
identically 1 but for elliptical orbit the ratio is [1 + ecos(A — Ap)]/ (1 — e^). This means, for 
an elliptical orbit the effective optical depth of the scattering region as seen by the source is 
increased by a factor of h{\ e) given in section 3 . For a circular orbit, e=0 and the effective 
optical depth and hence the degree of polarization becomes independent of the longitude. 
When A = Ap, i.e., when the planet is at the periastron, the effective optical depth as seen 
by the source is the maximum yielding maxiumum amount of polarization. 

As shown by Sengupta &: Krishan (2001); Sengupta (2003); Sengupta & Kwok (2005), 
the oblateness of the atmosphere plays a crucial role in determining the degree of polarization. 
The net degree of polarization is calculated by integrating the polarization at all points over 
the planetary disk. Since, we take the number density of scatterers distributed about some 
symmetry axis, the positive part of polarization gets cancelled out by the negative part . Now 
as oblateness increases, the net cancellation decreases owing to the departure from spherical 
symmetry. Consequently, the disc integrated polarization would vanish if the geometry is 
perfectly spherical and departure from sphericity would give rise to non-zero polarization. 
An increase in oblateness means more asymmetry and hence more net non-zero polarization. 
Figure 1 shows that if the oblateness of the planet is as small as 0.003 (similar to that of 
Earth) then the degree of polarization is of the order of 10~^ for the mean grain size of 
forsterite 10.0 jim. The degree of polarization estimated by Seager et al. (2000) for similar 
planets is of the order 10~^ for grain size 1-10 iim. This is because of the fact that multiple 
scattering reduces the amount of polarization by a few orders of magnitude as discussed in 
section 2. However, if the oblateness of the object is increased, the degree of polarization 
increases substantially as shown in figure 2 and could be at the range of present detectibility. 
In figure 3 we present the degree of polarization by forsterite grain of mean radius 10 iim 
for different values of the planetary oblateness. 

Figure 6 shows the polarization profiles of star with class II planet that orbits at a 
distance of 1-2 A.U and hence should have atmospheric temperature about 250 K. This 
type of planet should have a tropospheric cloud layer of water along with lesser amount of 
methane and ammonia. 

The amount of polarization is although dependent on the optical property of the con- 
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densates, the size of the condensates play the dominant role in determining the amount of 
polarization. According to the adopted chemical equilibrium model by Cooper et al. (2003), 
the particle number density increases with the decrease in grain size. As a result, the amount 
of polarization increases substantially when the mean particle size is decreased. Therefore, 
high degree of polarization might arise in a less oblate planet if the grain size is sufficiently 
small. On the other hand, large amount of polarization can be produced if the grain size 
is large but the oblateness is high as shown in figure 6. In figure 7, we present the degree 
of polarization caused by water particulates with different mean size. However, in order to 
understand the effect of grain size, we have kept the T-P profile the same as that calculated 
by Sudarsky et al. (2003) with mean water particulate size 5 iim. 

It may not be possible to infer from the polarization profile, the specific nature of the 
condensates, the oblateness of the planet or the orbital inclination angle as the amount of 
polarization would arise by a combination of all these parameters. However, if the distance 
of the planet from the primary star is known, then the nature of the condensate can be 
determined. Exact understanding on the chemical equilibrium may provide the size of the 
condensates. As a result, a combination of the oblateness and the inclination angle of the 
planet can be determined from the polarization data which in turn can provide information 
on the mass and the rotational velocity of the planet. However, the periastron angle and 
the orbital eccentricity can be inferred, irrespective of the amount of polarization, from the 
phase angles at which the maximum and the minimum amount of polarization arise. 



7. CONCLUSION 

The important message that is conveyed in this paper is that single scattering in a rota- 
tionally induced oblate planet may cause sufficient amount of linear polarization that could 
be detected by using the available instrumental facilities. The amount of linear polarization 
would be sufficiently large if the planet is highly oblate and the mean size of the condensates 
is not greater than a few micron. If detected, the polarization profile will provide informa- 
tion on the oblateness and orbital inclination of the planet and hence on the spin period 
and mass of the planet. Further, the time dependent polarization profile would provide the 
orbital period and the eccentricity. Our present discussion is confined to stars with single 
planet. The case for more than one planet orbiting a star at different inclination and phase 
angle is worth investigating as that should change the amount and the periodic nature of 
the variability of the polarization. The prsent investigation is aimed for the polarimetry of 
already known extra-solar planets. Since, the orbital separation and the emitted fiux from 
the parent star are known, one can employ the theoretical inferences about the composition 
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and location of the condensates in the planetary atmosphere. However, in general, the de- 
gree of polarization should depend on the evolution of the planets. Baraffe et al. (2003) have 
discussed on the time evolution of extra-solar planets and showed how the radius and hence 
the surface gravity as well as the effective temperature of an irradiated planet changes with 
time. Since the formation, chemical evolution and the geometrical location of condensates 
change with the change in surface gravity and atmospheric temperature profile, the degree 
of polarization, wheather by single scattering or by multiple scattering, would also alter with 
time. Therefore, time evolution of the polarization profile is worth investigating. 

We are thankful to the referee for several constructive comments and useful suggestions 
that have improved the clarity and the quality of the paper. This research has made use of 
NASA's Astrophysics Data System. 
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Fig. 1. — Degree of linear polarization at Bessel R band caused by forsterite grain with mean 
radius Tq = lOfim. The polarization is plotted as a function of the fractional orbital period 
t/P where t is time and P is the orbital period. Models for planet with oblateness f=0.003 
and circular orbit. From bottom to top at t/P = 0.25, the curves represent polarization 
with inclination angles i = 90°, 60°, 45", 30° and 0° respectively. 
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Fig. 3. — Degree of polarization at Bessel R band as a function of oblateness f. Circular 
orbit with inclination angle i = 90°. The linear polarization by planet with forsterite with 
ro = 10/xm is calculated at the fractional orbital period t/P = 0.0. 
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Fig. 4. — Same as figure 1 but for planets orbiting in elliptical orbit with eccentricity e = 
and the epoch of periastron passage Tq/P — 0. 
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Fig. 5. — Degree of polarization by planet with f=0.065 and Tq/P = 90" but different orbital 
eccentricity, t/ P is the fractional orbital period, t being time and P the orbital period. Solid 
lines from bottom to top at t/P = 0.5 represent models with e =0.0, 0.1, 0.2, 0.3, and 0.4 
respectively and with i = 90° while the dashed lines represent that with i = 30°. 
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Fig. 6. — Degree of polarization as a function of fractional orbital period t/P at Bessel R 
band due to water condensates, (a) Models for planet with f=0.003, e = 0.4 and mean 
particle radius Tq = 5.0 fim. Solid lines - models with Ap = 60° {Tq/P = 0.67); dashed lines 
- models with Ap = 120° {Tq/P = 0.42). Both the solid and the dashed lines from bottom 
to top at t/P — 0.25 represent polarization with i — 90°, 60°, 45°, 30° and 0° respectively, 
(b) Same as (a) but solid lines - models with tq = 1.0//m, f=0.001 and Tq/P — 0.25; dashed 
hnes - models with tq = 4.0//m, f=0.01 and Tq/P = 0.75. 
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Fig. 7. — Degree of polarization at Bessel R band as a function of mean grain diameter 
do- Circular orbit with inclination angle i = 90°. The linear polarization by planet with 
oblateness f=0.001 having water is calculated at t/P = 0.0. 



